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Abstract. We study the coupled Einstein-Yang-Mills-Dilaton (EYMD) equations for a Friedmann- 
Lemaitre universe with constant curvature k = 1. Our detailed analysis is restricted to the case 
where the dilaton potential and the cosmological constant vanish. Also assuming a static gauge field, 
we present analytical and numerical results on the behavior of solutions of the EYMD equations. 
For different values of the dilaton coupling constant we analyze the phase portrait for the time 
evolution of the dilaton field and give the behavior of the scale factor. It turns out that there are 
no inflationary stages in this model. 



1 Introduction 



For a long time the matter used in most cosmological models was described in terms of ideal 
or viscous hydrodynamics. This is certainly adequate at least back to the time of nucle- 
osynthesis, because — to a first approximation — we are then allowed to regard the matter as 
a mixture of pressure-free matter (baryons, WIMPS, etc..) and radiation (photons, neu- 
trinos). At very high energies the fluid description has to be replaced by field theoretical 
models (and, when approaching the Planck scale, by yet unknown theories of matter). 

*This work was supported by the Swiss National Science Foundation. 
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With the success of gauge field theories, and especially in connection with various pro- 
posed infiationary scenarios field theoretical models became very popular in cosmology. 
A description of topological defects, which could be created in cosmological phase transition, 
requires also field theoretical models which allow for spontaneous symmetry breaking. 

Inflationary stages appear more or less naturally in many models, but most of those 
which lead to successful inflation are very phenomenological in nature and are not based 
directly on "fundamental" fleld theories. Since inflation is a very attractive general idea 
and for other reasons as well, it is important to study the evolution of a large variety of 
cosmological models based on interesting fleld theoretical matter descriptions. In addition 
to Yang-Mills flelds, the existence of a dilaton fleld is theoretically well motivated on the 
basis of superstring theory and Kaluza-Klein theories [Q. 

In the present work we study the coupled Einstein- Yang-Mills-Dilaton (EYMD) equa- 
tions for a Friedmann-Lemaitre (Robertson- Walker) space-time. The basic equations for this 
model lead to a dynamical system of six flrst-order nonlinear ordinary differential equations, 
which we discuss with the modern techniques of the theory of dynamical systems. In view 
of the rather high dimension of phase space this goal is not easy to achieve. In a flrst step 
we therefore analyze in detail only a restricted class of solutions. We do so by setting the 
cosmological constant and the dilaton potential equal to zero, and also require that the Yang- 
Mills fleld is static, which can be done consistently. In this case the phase space becomes 
three-dimensional. We hope that our analysis of this reduced system is at least methodically 
of some interest. An extension of the discussion to the full system is planned. 

In a flrst step we compactify the phase space such that the dynamical system can be 
extended to the enlarged space and determine the critical points. With one exception, all 
critical points are either located at inflnity of the original phase space or in the unphysical 
region. We analyze the nature of these singular points and the qualitative behavior of the 
flow in their vicinity. In general the equilibrium points are hyperbolic. However, for special 
values of the dilaton coupling constant there exists a critical point in the asymptotic physical 
region which becomes non-hyperbolic. For this case we use the theory of normal forms and 
apply the Poincare-Dulac theorem to get rid of the nonresonant nonlinear terms up to third 
order. Analytical formula which are valid near the critical points are given in all cases. 

The asymptotic behavior provides one with the initial conditions for a numerical study 
of the phase portrait. Special attention is given to the value of the dilaton coupling constant 
which is obtained from superstring theory. Unfortunately, and not unexpectedly, we do not 
flnd any inflationary stages. 

The paper is organized as follows. In section 2 we derive the basic equations. Their 
analysis for the above-mentioned assumptions is performed in section 3. We analytically 
determine the asymptotic behavior of the solutions near the singular points. Section 4 is 
devoted to numerical results. The behavior of the scale factor as well as the phase portrait 
of the dilaton fleld are given. A short summary concludes the paper. 
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2 Basic Equations 



The metric of the Friedmann-Lemaitre universe with constant curvature /c = 1 is given by 

g = dt'^ - a\t') h = a\t) [dt" - h], (2.1) 

where h is the standard metric on and t is the conformal time related to cosmic (proper) 
time by dt' = a{t) dt. 

We work in the "Einstein conformal frame" in which the action for the EYMD system is 
given by 

^-i/(-|;-4 + >)^-^<*)-^^^)"- <^-^' 

where F = dA + AaA, k, denotes the dilatonic couphng constant and rj is the volume form. 
A and (p = (j){t) are the SU{2) Yang-Mills potential and the dilaton field, respectively. Both 
are homogeneous and isotropic. 

Since the group manifold of SU{2) is S^, there is a natural ansatz for a homogeneous 
isotropic SU{2) gauge field involving a single function of time 0: 

A = fit) e, (2.3) 

where G is the Maurer-Cartan form on regarded as SU{2). Note that t and / must be 
considered dimensionless while the scale factor a has the dimension of a length. The canonical 
metric h on 5*^ can also be expressed in terms of 0: If we decompose the sM(2)-valued form 
according to 

= E (^'r, , (2.4) 
where tj = —icj [aj are the Pauli matrices), then the 6^ are left invariant 1-forms and 

h = 6ij 6' O . (2.5) 

A simple way to verify the normalization is to use the following parameterization of SU (2) 
in terms of x = {x^,x^, x^, x^) G S^: 

— + ix^ x^ — ix^ J ' 
The Maurer-Cartan form is then U~^dU and the 6^ are given by 

6' = -hilTiU-'^dU] (z = 1,2,3). (2.7) 

We define an orthonormal basis of one-forms by 

■=adt, 9' ■.= ae' (i = 1,2,3). (2.8) 



The metric (2J) then reads 

g = g^J^^^e\ ^^, = diag(l, -1,-1,-1). (2.9) 
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Variation of the action with respect to the metric yields the Einstein field equations 

G^, = 87TG%, + Ag^,, (2.10) 



with 



-00 







4:71 


[2a2 






4:71 


[2a2 



+ V{<j>) + ^JP + 4{f-fr]e-'-^ 



and Toi = 0. The (00) component of Einstein's field equations is the constraint equation 

H' + l= G{h' + ^r(0) + + 4(/^ - me-'^''} + (2.11) 

with the Hubble parameter H = a/a, where an overdot means the derivative with respect 
to the conformal time t. For the spatial components we have 

- [2H + H^ + 1]= - 2aV(0) + \[f^ + 4{f - /)']e-2'"^} - Aal (2.12) 

The dilaton and the YM equation are obtained from variations of the action with respect to 
(p and /, repectively, 

+ 2i70 + a2^ + ^[f -4(/^-/)V'«<^ = 0, (2.13) 
/_24/ + 4/(2/-l)(/-l) = 0. (2.14) 

The system is fully determined by the independent equations ( |2.11D , ( ^.13| ) and ( [2.141 ). In 
fact, one can (easily) show that Eq. ( |2.12|) is a consequence of the others. This is a refiection 
of the Bianchi identities. 

From Eq. ( p.l4| ) we see that the static gauge field solutions are given by / = 0, 1/2, 1. 
Only the case / = 1/2 is of interest, because / = corresponds to a vanishing gauge field 
{A = 0) and / = 1 is a pure gauge {A = Q). 

A systematic study of the complete system is rather involved. We confine ourselves to 
the special case of a static gauge field / = 1/2 and choose in addition = and A = 0. 
The set of differential equations ( |2.11| )-( p.l4| ) then reduces to a three-dimensional dynamical 



system. In order to analyze this reduced system, we use the same procedure as in 0. We 
determine the singular points (see page 15), including those lying at infinity of the phase 
space, and then find analytically the asymptotic behavior for the solutions near these points. 



3 Asymptotic Behavior 

Setting / = l/2, V^ = 0, A = and introducing a "proper" time coordinate t according to 

,~ dt , , 

dt = ^-, 3.1 
a{t) ' ^ ' 
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we obtain the transformed system of basic equations 

- G[l<j)l + -e 



Hr 



-2K,<p-\ 



where H = a^/a and the subscript t denotes the derivative with respect to this variable. 



(3.2) 
(3.3) 



(3.4) 



Equations (3^) and (|3.3| ) form a three-dimensional dynamical system in the phase space 
0f , H. We now introduce the following dimensionless variables: 



t 



H 



V 
k 

X 

y 



Gt, 

3 i~ 
G 

H 



G 



Then Eqs. ( |3.2| ) and (|373|) become 



G 



and the constraint 



Vn 

Zn 



-yz + 



-2kx 



(3.5) 
(3.6) 

(3.7) 



takes the form + \e 



a} with z = dfj/a. 



For the special case k = 0, these equations are equivalent to a one-dimensional mechanical 
system for a{ri) with the Lagrange function L = T — V = a^ — {a'^ — a^/i). Indeed, Eq. ( p.6|) 
immediately implies the relation y = k/a, where k is an arbitrary constant. From the 
constraint equation we obtain the corresponding conserved "energy" E = T + V = dl + 



a^/4) = /c^ > 0. The solution with respect to the conformal time is 



a(t) 
z{t) 

x{t) 



72 



1 + Vl + 64A;2 sin (2t - A) 



VI + 64P cos (2t - A) 
4 (H- Vl + 64A;2 sin(2t- A)) 



2 



In 



tan (t - A/2) + VTTeiF - 8A;' 
tan {t - A/2) + VTTeiF + 8A; 



+ const. t G (0, 7r/2 + A), 



where A is defined by 

X ■ ( ^ 

A := arcsm 



VI + 64A; 
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If K 7^ 0, we use instead of x the function 



« = ^ > 0. (3.8) 

Our system then takes the final form 

i/r, = -yz + Ku^, (3.9) 

2 2 r, 2 

= z — u — 2y , 

with the constraint 

+ ?/2 - = (3.10) 
where we recall that z is the Hubble parameter with respect to rj, 

z=^. (3.11) 
a 

The system of equations (|3.9|) - (|3.11|) possesses the following symmetries: 

a) u ^ —u h) rj ^ —rj c) a ^ —a d) y ^ —y 
y -y -K 

z —z 

According to a) the plane m = is a two-dimensional invariant subspace of the phase space; 
it separates solutions with u > from those with m < 0. Trajectories can reach this plane 
asymptotically, but they never penetrate it. We are not interested in solutions with m = 0; 
for these the change of variables given by Eq. (^.8|) would no longer be valid. 



Equations (|3.8|) and ( |3.10 ) imply that the trajectories of our model lie in the region of 
phase space exterior to the cone + — = where m > (see Fig. 1). 

We also note that the sections of the trajectories lying in the upper half of the phase 
space {z > 0) correspond to expansion {H > 0) and those in the lower {z < 0) to contraction 
{H < 0). When a trajectory intersects the plane z = {H = 0), the scale factor a attains 
its maximum, because Zri = —{v^ + 2y^) is strictly negative. 



3.1 The Singular Point at the Origin 

Looking for equilibrium states of the dynamical system (p.9|) we immediately see that the 
origin {u,y,z) = (0,0,0) is the only singular point not lying at infinity. However, this 
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Figure 1: Physical domain. 

critical point is not of much interest, because there are no trajectories originating from it 
which extend to the exterior of the cone nor does any trajectory (coming from the exterior 
of the cone) ever reach the origin: Suppose there would be a trajectory emerging from the 
origin {u,y,z) = (0,0,0). Along this, z = a^/a has to be negative, since Zr, = —{y"^ + a?) 
is strictly negative. It follows that if a,, > (a,, < 0) along the trajectory, then d must be 
negative (positive). This is a contradiction because d vanishes at the origin. With similar 
arguments one can prove that no trajectory can reach the origin from the exterior of the 
cone. 

The repelling behavior of the singular point {u,y,z) = (0,0,0) is nicely illustrated by 
considering the case where k = 1. Then the plane z = —y is an invariant subspace and the 
phase portrait of the trajectories lying on this plane is given in Fig. 2. Furthermore, we can 
give the explicit solution. On the invariant plane z = —y our system of Eqs. (|3.9|) reduces to 



Uri = UZ, (3.12) 
Zri = —z^ — U^, (3.13) 



and the constraint equation becomes = a^. Eq. ( p.l3|) is equivalent to the equation of 



motion of a one-dimensional mechanical system with the Lagrange function L = T — V = 
— The corresponding conserved "energy" is £^ = T + = a^/2 + > 0. The 

solution with respect to the conformal time is 

a(t) = u{t) = {2 Ve sin {V2t))^, 
\fE cos 

x{t) = -iln [s/E sin(y2t)] te{0,7i/V2). 

In the limit t ^ 0, the scale factor a vanishes and the equation of state tends to e = p (stiff 
matter), where e = Too and p = ^{Tu + T22 + T33). At t = -\/2f , d attains its maximum and 
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the equation of state becomes e = 3p which corresponds, for ordinary matter, to a extremely 
relativistic gas. 




-300 -200 -100 100 200 300 

y 



Figure 2: Phase portrait of the invariant plane z = —y. 



3.2 The Singular Points at Infinity 

In order to find and investigate the singular points at infinity it is convenient to compactify 
the phase space and complete it by the infinitely distant boundary u"^ + y"^ + z"^ = oo. We 
do so by going over from the Cartesian coordinates u, y, z to the spherical coordinates 

u = r cosipsini), y = r sincp sini), z = rcos'd, (3-14) 

and by performing a subsequent transformation of the radius according to 

^ (0 < p < 1). (3.15) 



1-p 

Also introducing a new time r defined by 

dr]={l-p)dT, (3.16) 

and using the variables p, 'd, (p, r, our system ( p. 91 ) assumes the form 

Pr = p^(l — p) cos'i9[cos(2'i9) — 2 sin^ sin^ yj], 
= — p sin 

79cos(279)[l + sinV], (3.17) 
ifr = p COS (p[k sin i} — simp COS i)]. 



For < p < 1, Eqs. ( pl9| ) are equivalent to the system ( p.l7| ), but the latter admits a 



smooth continuation to the boundary p = 1. Thus, in the variables 0<p<l,0<'i9<7r, 
< (fi < 271, the phase space of the system ( |3.17| ) is compact. Each point of the original 



phase space is mapped along its radius vector to a certain point inside the unit ball, while 
the points + y"^ + z"^ = oo a.t infinity are mapped to its surface p = 1. From Eqs. ( |3.17| ) it 
is readily seen that on this surface the system has several singular points, their number and 
nature depending on the value of k. Due to the symmetry transformation d) it is sufficient 
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P' 



Figure 3: Singular points at infinity. 

to consider only positive values of k (i.e. k > 0). There are at least six singular points which 
will be denoted by (see Fig. 3) 

P(^ = 0), K^m.ip^) = (7r/4,37r/2)), K2m,V2) = (V4, V2)), 
P'(^' = 7r), ir;((^;,y,;) = (3V4,V2)), ir^((^'2,^2) = (3V4,3V2)). 

For fi: < 1 we have the four additional points 

K2.{{^^,^^) = (7r/4,arcsinK)), -ft'4((^?4, = (7r/4, vr - arcsin k)), 
K'sii^s,^^) = (37r/4, -arcsin^)), K'^{{^'^, ^'^) = (37r/4, vr + arcsin^)). 

In the limit k — 1, these coincide with K2 and K2 respectively. The symmetries a) and b) 
imply that there are only four essentially different points: P, Ki, K2 and K^. Since the 
singular point P at the north pole lies in the interior part of the cone, it is not of physical 
interest. 

Before discussing the remaining singular points we would like to make the following 
remark: For the three points Ki {i = 1,2,3) there exists a separatrix which goes from Ki 
straight down to the origin (0,0,0), moving on the surface of the cone. In fact, Eqs. ( ^.17| ) 
have admit solutions of the form {{}{t),(p{t)) = {{}i,ipi), where {{}i,ipi) denotes the angular 
coordinates of Ki, along which we have Pr < 0. This is not contradictory to the statements 
made in section 3.1, because the scale factor a vanishes along these separatrices. 

3.2.1 The Singular Point Ki 

In order to study the asymptotic behavior of the solutions near Ki we first perform a shift, of 
the variable, defined by 5p = p — 1, 6i} = i) — -^i and 6(p = ip — (pi. In these new coordinates. 
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the linearization of Eqs. ( |3.17| ) at Ki gives 





( K 








z): 












^ 









(3.18) 



with the strictly positive eigenvalues \p = \^ = 2\/2 and = V2/2 (1 + k). Hence, 

the singular point Ki is hyperbolic p. 69] and the leading terms near Ki are given by 

(MM = -{Sp)oe'^\ (5^)(r) = (5^),e^-^ 

{5ip){r) = (Moe^-^ (r^-oo), (3.19) 

with arbitrary integration constants (5p)o > 0, {S{})o > and {Sip)o > 0. Integration 
of Eq. (|3.16|) leads to the relation e^'''^ = (^^^; i-e., emergence from Ki corresponds to 
increasing rj, starting from 7] = 0. 

For the Cartesian coordinates {u, y, z) we obtain (indicating only the leading terms as 
V 0+) 



u{vi) = {6^)c 



(Mo. 



y{v) = --, (3.20) 

V 

1 



z{ri) 



V 



From Eq. (|3.11| ) it is readily seen that near Ki the scale factor follows the law d{r]) = C rj, 
where C > is a constant, fixed by the constraint ( |3.10|) . 

Hence, to the point Ki there corresponds an initial cosmological singularity at a certain 
finite time which can always be chosen at t = 0. Near Ki we have 



where to > is an arbitrary constant. Emergence from Ki corresponds to increasing t, 
beginning from t = 0. The effective equation of state near Ki tends to e = p (i.e., stiff 
matter). 



3.2.2 The Singular Point K2 

In order to study the asymptotic behavior of the solutions near K2 we proceed in a similar 
manner. Linearization at K2 again leads to Eq. ( ^181) with the same eigenvalues Ap and 
A^, but with the changed eigenvalue A^ = \/2/2 (1 — k). We suspect that in this case the 
behavior may depend on the value of k. For k < 1 all eigenvalues are again strictly positive. 
Hence, the singular point K2 is hyperbolic and the leading terms in the expansion around 
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K2 are given by Eq. ( |3.19| ) with arbitrary constants {5p)o > 0, {6{})o > and {Sip)o < 



0. Integration of Eq. (|3.16|) leads to the relation 



T], i.e. emergence from K2 



corresponds to increasing 77, beginning from 77 = 0. 

For the Cartesian coordinates {u, y, z) we obtain in leading order as 77 — 0"' 

A. 



u{ri) 

y{v) 

z{r]) 



-(Mo 



(Mo 



V 



1 

V 

1 



(3.22) 



Near K2 the scale factor follows the law d{ri) = C r], where C > is a constant fixed by 
the constraint (|3.10|) . Hence, if k < 1, K2 corresponds to an initial cosmological singularity 
at a finite time t = 0, say. Near K2 we have 



H 



1 

2t' 



3 1 
4Gt' 




4G 4,' 



(3.23) 



where to > is an arbitrary constant. Emergence from K2 corresponds to increasing t, 
starting from t = 0. The equation of state tends to e = p (i.e., stiff matter). 

If K > 1, the singular point K2 is also hyperbolic, however, X^p is strictly negative. The 
outgoing trajectories necessarily lie on the invariant plane u = (i.e: {S(p)o = 0) and are 
thus not of interest. 

The case k = 1 is more involved because the singular point K2 is non-hyperbolic. Terms 
of higher order become important. Up to third order the expansion in a sufficiently small 
neighborhood around K2 is 



(M. 



^6p + V26p'+^-^6p6^ + y^6p' 



2V2 5d + 2^2 5^"^ + 2^2 5 p6d 
7^/2 



(3.24) 



a/2 



We can now use the Poincare-Dulac theorem 
in Eq. 



|, p. 191] to classify the nonlinear terms 
^.24| ) into resonant and non-resonant ones. For the non-resonant terms of order n 
there is a polynomial change of coordinates of degree n so that they are transformed into 
terms of at least order n + 1. This is not the case for the resonant terms. The polynomial 
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change of coordinates is found by solving the so-called homological equation. Performing 
the polynomial change of coordinates on Eq. ( p.24|) , the non-resonant terms of degree 2 and 



3 become of higher order and are thus neglected. Retaining only terms up to third order, 
we obtain the set of equations 



dvp V2 V2 ^ 



dr 



v^v , (3.25) 
dv^ _ _\/2 3 

where {6p, 6'^, 6ip) = {vp, v^, v^) + polynomes of degree 2 or higher in {vp, v^, v^). The minus 
sign in the last equation implies that for k = 1 there are no trajectories which emanate from 
K2 into the exterior of the cone. Indeed, Eqs. ( p.25| ) can be solved analytically. We obtain 

= (^^^"''"'^ • 

t \ ^'&o 4/ N 2V2T 

where Vpo,v^o and t>^o 7^ are arbitrary constants. The + (— ) sign corresponds to the case 
Vipo > {v^po < 0). Hence, for non- vanishing constants Vpo and f^o? we have 

v^{r) ^ 0, |fp(r)|^oo, \v^{r)\ ^ 00 (r ^ 00). 



3.2.3 The Singular Point Kg 



The singular point K3 only appears for n < 1. Linearization at leads to 





/ 6p 

" ' — (1 + K^) and /i^ = 
i^s is thus hyperbolic and the leading terms near it are given by 



/ip U 

/i^ 

\pl'k\J\ — /i^ 




with the eigenvalues /Xp = \/2/2k'-^ jj,^ 



(3.26) 

v^/2(l - k^). The point 



(M(r) 
(M(r) 



{6^)06^'^^ (r 



-00 



(3.27) 



with arbitrary constants (5p)o > and {S'd)o > 0. Note that the general solution of 
would be given by 

(MM 
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However, since /x^ is negative the constant {5ip)o must vanish. 



Integration of Eq. (|3.16| ) leads to the relation e^'''' 



j-^^^ ?7, i.e., emergence from 
corresponds to increasing r/, starting from r/ = 0. In terms of the Cartesian coordinates 
(m, I/, z) we obtain 



u{ri) 



T] * 



3 - -2 



K 



3 + K 



Z-2 



mo 



-, M 



Z{7]) 



^«^ + «^ 't— ^ 5^ o 

?7 2 S + 

V 2 [(5p), 



.(Mo, 



(3.28) 



-, fill 
p-p p±-i 
fjPp ^ 



with ^ 

/ip 



2 1 + K 



> 4. From Eq. (|3.11D it is readily seen that near the scale factor 



behaves like 0(77) = C i]^^'^ , where C > is a constant, fixed by the constraint (p.lO|) . 



The point therefore corresponds to an initial cosmological singularity at a finite time 
t = 0, say. Near we have 



H = {l + k 



K 



3 1 



1 + ^2 V G t ' 



1+^2/1 14), 



(3.29) 



where to > is an arbitrary constant. Emergence from K^, corresponds to increasing t, 
starting from t = 0. 



4 Numerical Results 



The asymptotic behavior of the solutions near the singular points provides one with the initial 
conditions for solving the differential equations numerically. We will recover the repelling 
{Ki,K2^K^) and attracting {K[, K'^^ K'.^) characteristics of the singular points, as well as 
their dependence on the value of R. Special emphasis is given to the value R = ^/3 which is 
predicted by a model obtained from superstring theory 
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4.1 Emergence from K] 



The asymptotic behavior (|3.2CI|) of the solutions near Ki provides proper initial conditions 
at a given time t] = rjo- The constraint is taken into account by fixing the constant C at the 

time rjo- We can now produce plots for different values of R and k := {Sip)o 



l+k 
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k = 1/2 : Fig. 4 shows the behavior of the scale factor a for different values of the constant 
k. The phase portrait of the dilaton field is given in Fig. 5. Since k <1, each singular point 
K[, K2 or K'^ is attractive. From Fig. 6 it can be seen that the plots actually correspond to 



400 




400 




0.015 



Figure 4: The scale factor with respect to the time r] and the conformal time t respectively. 




Figure 5: Phase portrait of the dilaton field. 



37r/4 




-7r/2 + e 


= 5/ \ 












= 1 


-7r/2 









0.015 



0.015 



Figure 6: The behavior of d and (p with respect to 77. The value of £ is £ = 8 • 10 ^. The 
trajectories are therefore close to the plane u — Q. 



trajectories which end at i^Tg. We have also computed numerically •& and </? as a function of 
time for all subsequent cases. This enable us to identify the points at which the trajectories 
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end. These plots can be found in the Diploma thesis of one of us 
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K = : Fig. 7 shows the behavior of the scale factor a for different values of the constant 
k. The phase portrait of the dilaton field is given in Fig. 8. Since for k > 1 only K[ is an 
attractive point, we expect that the plots correspond to trajectories which end at K[. This 
is confirmed by our numerical computation. 




Figure 7: The scale factor with respect to the time t] and the conformal time t respectively. 




X 



Figure 8: Phase portrait of the dilaton field. 



4.2 Emergence from K2 

We proceed in a similar manner. Emergence from K2 only appears ii k < 1. The constant 
k is defined by k := -{6(^)0 (see Eq. (^)). 

k = 1/2 : Fig. 9 shows the behavior of the scale factor a for different values of the constant 
k. The phase portrait of the dilaton field is given in Fig. 10. The plots correspond to 
trajectories which end at the singular point K[. 
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Figure 9: The scale factor with respect to the time 1] and the conformal time t respectively. 
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Figure 10: Phase portrait of the dilaton field. 
4.3 Emergence from K3 

The singular point only exists if k < 1. We define the constant k hj k := [5^) 
(seeEq. (pg])). 




Figure 11: The scale factor with respect to the conformal time t. 
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k = 1/2 : Fig. 11 shows the behavior of the scale factor a for different values of the constant 
k. The phase portrait of the dilaton field is given in Fig. 12. The plots correspond to 
trajectories which end at K[. 
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Figure 12: Phase portrait of the dilaton field. 



5 Concluding Remarks 

For the special case of a static gauge field and V = A = we have performed a complete 
analysis of the coupled EYMD equations for a Friedmann-Lemaitre universe with constant 
curvature k = 1 (closed model). For this situation, we did not recover any inflationary 
stages. It might therefore be interesting to extend our discussion to the more general case 
of a dynamical gauge field, perhaps with a non- vanishing cosmological constant A and/or 
a dilaton potential V . This is, however, a considerably more involved task, because we 
encounter then six coupled first-order nonlinear differential equations. 

The analysis of the coupled EYMD equations for a Friedmann-Lemaitre universe with 
constant curvature A; = (fiat model) and k = —1 (open model) are also of interest, since 
related considerations [0, |ll[] show that inflationary stages occur for trajectories which come 
close enough to one of the separatrices found for k = 0. 
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